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Using  both  sea  water  data  and  research  tank  water  data,  cou¬ 
pled  with  the  theoretical  approach  of  Wells,  calculations  of  the 
modulation  transfer  function  (MTF)  for  the  water  medium  are 
performed  via  use  of  the  Wells  Decay  Function,  Various  approxi¬ 
mate  relationships  are  used  to  derive  the  value  of  h,  the  total 
scattering  coefficient,  which  enters  into  the  MTF  calculation. 
The  approximated  b  values  are  used  to  compute  the  water  MTF. 
These  estimated  MTF  calculations  are  compared  with  experi¬ 
mental  data  MTF  calculations  for  which  the  actual  water  value 
of  b  was  derived.  Error  assessments  and  recommendations  for 
the  various  b  relationships  used  in  the  MTF  computation  are 
made. 
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Introduction 

The  effect  of  suspended  particles  in  a  water  medium  on 
the  transport  of  light  is  normally  dichotomized  into  two 
processes — namely,  the  dual  effects  of  absorption  and 
scattering.  Absorption  acts  simply  to  remove  photons 
that  might  have  contributed  to  image  formation.  These 
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absorption  attrition  effects  might  be  overcome  for  a  given 
system  by  increasing  the  source  flux  or  by  using  a  more 
sensitive  receiver.  Scattering,  however,  in  contradistinc¬ 
tion,  is  more  pernicious  through  two  allied  processes. 
Scattering  can  produce  a  foreground  veiling  glare  that 
reduces  the  effective  contrast  of  the  target  with  respect  to 
the  background,  and,  by  causing  a  redistribution  of  the 
trajectories  of  image-forming  photons,  scattering  scram¬ 
bles  the  information  content  originally  presented  by  the 
target  and  background. 

The  modulation  transfer  function  (MTF)  is  used  to 
characterize  the  effects  the  medium  has  on  the  passage  of 
image-forming  light.  It  is  not  surprising  that  values  of 
optical  parameters  that  describe  the  scattering  effects  of 
the  water  medium  are  significant  variables  in  the  MTF 
expression.  Wells1,2  has  derived  a  transformation  that 
converts  the  medium  volume  scattering  function  (VSF), 
/3(W)  [  /3(H)  represents  the  fraction  of  radiant  power  in  a 
beam  that  is  scattered  by  angle  (-)  per  steradian  per  unit 
length  along  the  beam],  into  a  medium  MTF  expression 
via  the  use  of  a  'decay'  function.  The  MTF  expression  so 
derived  reouires  the  value  of  the  total  (volume)  scatter¬ 
ing  coefficient,  b.  The  value  of  b  is  related  to  the  VSF  by 
the  functionality 

b  =  0(f)) sin  OcfO,  (1) 

129 


stated  here  in  spherical  coordinates  and  assuming  az¬ 
imuthal  symmetry.  Figure  1  provides  a  diagram  showing 
the  scattering  process  and  identifying  the  scattering  an¬ 
gle  0.  Azimuthal  symmetry  of  the  scattering  process  is 
assumed  for  scattering  in  water  media. 

Normally,  6  is  not  measured  directly,  because  of  either 
the  lack  of  necessary  instrumentation  or  simply  the 
time-consuming  nature  of  measuring  (3(0)  for  as  many 
angles  as  practical  while  holding  station  at  sea.  Rather, 
the  value  of  b  is  computed  by  using  other  more  easily 
measured  optical  coefficients.  Preisendorfer,3  Honey,4 
Wilson,5  Timofeyeva,6  and  Morel7  provide  certain  spe¬ 
cific  relationships  that  can  be  used  to  determine  b  for  a 
given  water  mass.  These  various  relationships  are  based 
on  sundry  data  set  interpolations  or  fitting  procedures. 
The  data  sets  and  the  conditions  of  their  collection  and 
subsequent  processing  provide  the  differences  seen  in  the 
various  expressions  for  estimating  b.  It  would  be  useful 
to  assess  the  value  of  these  various  expressions  for  the 
problems  associated  with  imaging  through  ocean  water. 

The  objective  of  this  paper  is  to  utilize  Wells’  results  to 
compute  the  medium  decay  function,  Z)(D,  for  a  homoge¬ 
neous  water  mass  (from  which  is  easily  determined  the 
germane  MTF  for  a  stated  range,  R),  using  both  the 
value  of  b  computed  with  the  above-mentioned  relation¬ 
ships  and  experimentally  derived  values  of  b.  The  result¬ 
ing  decay  function  expressions  will  be  compared  to  pro¬ 
vide  a  measure  of  the  errors  introduced  by  using  the 
various  estimation  methods.  Finally,  comments  will  be 
made  on  which  model  appears  to  provide  the  most  accu¬ 
rate  estimation  of  6  for  different  water  mass  types. 

Background  and  Methods 

The  medium  MTF  is  the  sine  wave  frequency  ampli¬ 
tude  response  of  the  water  path.  Wells’  transformation12 
of  the  VSF  into  a  medium  MTF  assumes  the  small  angle 
approximation  and  assumes  that  the  suspended  particles 
are  larger  than  the  optical  wavelength. 

Wells’  expression  for  the  medium  MTF  can  be  written 
in  the  form2 

T(T,R)  =exp[-Z?(r)R],  (2) 

where  R  is  the  range  of  interest  and  £>(D  is  the  ‘decay’ 
function,  with  F  the  spatial  frequency  in  cycles/radian. 
The  decay  function  for  non-absorbing  waters  may  be 
written  as2 

D(r)  =  6  -  <?(['),  (3) 

where  Q(D  is  the  frequency  space  expression  of  (3 (0),  or 
as  the  Hankel  transform  of  the  VSF.  Q(  D  provides  for  a 
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Figure  1.  Scattering  of  an  incident  light  beam  by  an  elementary 
volume  element. 
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recovery  of  power  at  lower  spatial  frequencies,  where  the 
scattering  is  too  small  to  create  blur  for  target  con¬ 
stituents  whose  relative  dimensions  are  large.1  Here  it  is 
assumed  that  the  scattering  function  is  highly  peaked  in 
the  forward  direction.  Thus,  the  small  angle  approxima¬ 
tion  may  be  used.  In  the  oceans,  this  is  a  very  robust 
assumption.  The  particles  present  are  generally  much 
larger  than  the  wavelength  of  light  incident  upon  them, 
and  their  refractive  index  approaches  that  of  the  water 
itself  Hence,  the  refraction  that  occurs  is  small,  and, 
generally,  diffraction  dominates  the  scattering  process. 

It  is  important  to  notice  that  the  decay  function  is 
independent  of  the  range,  R.  For  collimated  light,  Q(D  is 
expressed  by 1 

Q(  T)  =  2tt["  f}(0)Jn(a)0  d6,  (4) 

Jo 

where  (3(  0 )  is  the  VSF  and  J0(a)  is  a  zeroth  order 
Bessel  function  with  a  =  27rr@. 

It  can  be  shown,1  for  nonabsorbing  waters,  that  Q{  D 
monotonically  decreases  as  spatial  frequency  increases 
(i.e.,  as  the  potential  for  image  detail  increases),  and,  in 
the  limit  of  high  spatial  frequency,  Q(  F )  approaches  zero. 
In  this  region,  /XT)  approaches  6.  Hence,  the  diminution 
of  fineness  of  resolution  of  the  image-forming  light  over  a 
range,  R,  is  due  to  an  exponential-type  decay  of  the 
relevant  spatial  frequencies,  which  is  controlled  by  the 
value  of  b.  If  absorption  effects  are  included,  then  the 
effects  of  absorption  must  also  be  taken  into  account,  so 
that  D(D  may  be  written  in  the  form 

D(D  =  a  +  6  -  Q(r),  (5) 

where  a  is  the  absorption  coefficient. 

Table  I  provides  the  relationships  used  to  estimate  b 
knowing  other  optical  parameter  values  for  the  water 
mass  in  question.  In  Table  I,  c  represents  the  beam 
attenuation  coefficient,  k  the  diffuse  attenuation  coeffi¬ 
cient,  and  Chi  the  concentration  of  chlorophyll-a. 

The  scattering  data  used  in  the  computations  were 
collected  in  a  large  research  tank  facility  during  laser 
backscatter  experiments  conducted  at  the  Naval  Re¬ 
search  Laboratory”  (NRL).  The  tank  water  was  cleaned 
by  continuous  pumping,  filtering,  and  sedimentation. 
Then,  beginning  from  the  initial  clear-water  state,  the 
water  was  made  increasingly  more  turbid  in  a  controlled 
sequence.”  The  water  was  closed-circulated  through  the 
tank  to  keep  the  turbidity  agent  entrained.  This  system 
insured  that  the  water  in  the  tank  could  be  assumed 
homogeneous.  Ancillary  supporting  optical  measure¬ 
ments  were  made.  These  measurements  included  the 
beam  attenuation  coefficient,  c  (which  is  the  sum  of  the 
absorption  coefficient,  a,  and  the  total  scattering  coeffi¬ 
cient,  b),  the  diffuse  attenuation  coefficient,  k,  and  the 


TABLE  I.  Relationships  Used  to  Compute  b* 


Name 

Relationship 

Parameters 

Honey 

b  =  6/5  <  (c  k ) 

C.  k 

Wilson 

6  =  1.18v  (c  -  k) 

c.  k 

Preisendorfer 

b  =  c  3/4  x  k 

c.  k 

Timofeyeva 

klc  =  [0  19  x  (1  iv)]"  ' 

c,  k{w  =  b  i  c) 

Morel 

b  =  0.30  x  <Chl)°  62 

Chi 

*c  represents  the  beam  attenuation  coefficient,  k  the  diffuse  attenuation 
coefficient,  and  Cht  the  concentration  of  chlorophyll 
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Figure  2.  Measured  VSF  from  sample  experimental  data. 


0  1  0  20  30  40  50  60  70  80  90  1  00 

Spatial  Frequency  (cycles/rad) 


Figure  4.  Wells’  decay  function  for  an  intermediate  value  of  b 
from  test  tank  data. 


VSF,  0(0).  Figure  2  provides  an  example  of  an  experi¬ 
mental  VSF  measured  during  the  experimental  runs. 

The  VSF  measurements  were  fitted  by  a  modified 
Heyney-Greenstein  function9  so  that  values  of  the  total 
scattering  coefficient,  6,  could  be  computed,  using  Eq.  (1) 
for  the  different  water  mass  turbidities.8  The  relevant 
wavelength  for  the  tank  data  used  in  this  work  is 
528  nm. 

Results  and  Discussion 

As  shown  in  Eq.  (2),  the  MTF  can  be  determined  by  the 
negative  exponential  of  the  product  of  the  decay  function 
and  the  range.  Hence,  to  make  relevant  calculations 
independent  of  the  range,  the  decay  function  will  be 
computed  for  the  germane  b  values.  The  MTF  can  then 
be  found  easily  by  use  of  Eq.  (2). 

Figures  3,  4,  and  5  provide  the  decay  function  for 
different  experimental  and  estimated  values  of  b.  Here  it 
is  understood  that  the  values  of  the  decay  function  have 
not  been  normalized  by  the  argument  value  of  the  Bessel 
function.  Moreover,  to  provide  for  absolute  attenuation,2 
the  decay  function  now  includes  the  effects  of  absorption. 
Figures  3,  4,  and  5  provide  six  curves  each.  The  solid  line 
is  the  decay  function  computed  with  the  experimentally 
derived  b.  The  other  curves  shown  provide  the  decay 
curves  for  b  estimated  using  the  relationships  given  in 
the  legends.  Note  in  these  plots  that,  in  the  limit  of  high 
spatial  frequency,  the  value  of  the  decay  curves  tends 
toward  the  sum  of  the  absorption  and  total  scattering 
coefficients,  a  +  b,  as  they  should.  In  the  limit  of  zero 
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Figure  3.  Wells’  decay  function  as  a  function  of  spatial  fre¬ 
quency  and  value  of  b  used. 


spatial  frequency,  the  decay  function  approaches  the 
value  of  the  absorption  coefficient.  This  occurs  because 
the  target  in  this  limit  is  relatively  large,  so  that  any 
scattering  would  not  create  blur  and  the  only  loss  would 
come  through  absorption  effects  alone. 

The  mean  bias  error  is  given  in  Table  II.  The  mean 
bias  error  was  computed  by  taking  the  mean  of  the 
difference  between  the  different  fc-estimator-value  decay 
functions  and  the  experimentally  derived  b  decay 
function. 

Essentially,  the  decay  curves  shown  signify  how  accu¬ 
rate  the  methods  for  estimating  b  are  when  compared 
with  actual  b  values.  Figure  6  provides  a  plot  of  b  versus 
c  values  from  the  research  tank  data.  Here  it  is  seen  that 
the  Honey/Wilson  (these  two  are  very  close)  relationship 
over  much  of  the  range  covered  by  the  plot  provides  the 
most  accurate  estimate  of  b. 

However,  the  values  of  b  for  the  tank  were  relatively 
high.  Perhaps  in  real  world  waters  with  lower  values  of 


TABLE  II.  Mean  Bias  Error  of  the  b  Relationships  Used 
Compared  with  Experimentally  Derived  b  Values 


Method 

Mean  Bias  Error  (r; ' 

Timofeyeva 

41.7 

Honey 

7.4 

Preisendorfer 

8  6 

Wilson 

69 

Morel 

40  1 

Computation  of  the  Medium  MTF 
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Figure  6.  Research  tank  results  for  estimated  and  experimental 
b  values. 


different  set  of  optical  conditions — much  like  a  new  lens 
element  of  an  optical  system.  Thus,  for  these  cases,  an 
MTF  may  be  computed  for  each  layer  and  the  results 
cascaded  to  provide  an  overall  medium  MTF.  The  b 
values  to  be  used  for  the  MTF  of  such  a  medium  must 
then  be  profiles  of  the  b  value,  that  is,  values  of  b  as  a 
function  of  depth.  Using  the  Honey/Wilson  or 
Preisendorfer  approach,  profiles  of  k  and  c  are  needed. 
These  may  be  obtained  with  standard  optical  instrumen¬ 
tation.  However,  the  Morel  relationship  requires  a  knowl¬ 
edge  of  the  chlorophyll  concentration  profile.  Often  this  is 
obtained  with  a  fluorometer  that  is  calibrated  to  read  the 
chlorophyll  concentration  values  in  the  water.  Kitchen 
and  Zaneveld10  have  recently  shown  that  the  chlorophyll 
maxima  do  not  necessarily  coincide  with  all  potential 
scattering  present,  because  of  the  possibility  of  seston 
(the  total  concentration  of  suspended  material  over  the 
depth  of  interest)  layer  scattering  at  other  depths.  Thus, 
if  the  Morel  relationship  is  used,  some  means  of  deter¬ 
mining  additional  scattering  layers  independent  of 
chlorophyll  concentration  appears  to  be  needed. 
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Figure  7.  Sargasso  Sea  measured  b  data  versus  estimated  b 
values. 


scattering,  a  different  method  of  estimating  the  value  of 
b  would  be  more  successful.  To  check  this  proposition,  an 
NRL  data  set  for  the  Sargasso  Sea  region  was  used,  for 
which  the  basic  quartet  of  irradiances,  upwelling 
scalar/cosine  and  downwelling  scalar /cosine,  were  mea¬ 
sured  at  a  wavelength  of  520  nm.  Figure  7  exhibits  the 
results  of  the  calculation  of  b.  It  is  seen  that,  in  these 
very  clear,  oligotrophic  waters,  over  a  small  range  of  c 
and  b  values  determined  by  the  data  set,  the  Preisendor¬ 
fer  relationship  appears  to  provide  the  most  accurate 
estimate  of  b  overall.  This  is  closely  followed  by  the 
Wilson  relationship. 

In  our  computations,  it  is  assumed  that  the  water  path 
of  interest  is  horizontal  (that  is,  a  plane-parallel,  strati¬ 
fied  medium  is  assumed).  If  slant  or  vertical  paths  are 
considered,  then  each  layer  of  the  medium  represents  a 


Summary 

Decay  functions  were  calculated  for  both  experimen¬ 
tally  determined  and  estimated  values  of  b.  The  relevant 
MTF  may  be  calculated  directly  from  the  decay  function 
curve,  given  the  range,  R,  of  interest  [see  Eq.  (2)].  A  table 
listing  of  average  error  over  the  b  values  used  was  given. 

Of  the  relationships  tried  (Table  I),  the  Honey/Wilson 
relationship  proved  to  give  the  most  accurate  estimation 
of  b  for  the  relatively  turbid  research  tank  data.  How¬ 
ever,  for  the  clear  waters  represented  by  the  Sargasso 
Sea  data  set,  the  relationship  of  Preisendorfer  gave  bet¬ 
ter  estimates  of  b.  A 
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